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When the flavour content of QCD is increased sufficiently, the theory develops a non-trivial infra 
red fixed point. Thus, for a number of flavours above a certain critical value, but not yet so 
high that asymptotic freedom is lost, QCD becomes a conformal field theory. The location of 
the lower limit of this conformal window has not yet been unequivocally determined. Using an 
improved lattice action, and exploiting modern algorithms allowing for larger lattices and lower 
quark masses, we have shown that the theory of QCD with eight flavours breaks chiral symmetry 
in the continuum. We present proof that the accompanying transition is thermal in nature and 
as a consequence, the conformal window of QCD can only start afterwards, corroborating recent 
analytical studies at the expense of older results. 
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1. The conformal window 

In seminal papers, Caswell [1] and subsequently Banks and Zaks [2] showed the existence 
of an additional zero of the beta function for non-abelian gauge theories with fermionic content 
in the fundamental representation and for certain values of the number of flavours Nf. Their re- 
sult suggests that QCD Uke theories with a number of massless flavours within a certain range are 
both asymptotically free and chirally symmetric at zero temperature. Long range dynamics will 
therefore exhibit the qualities of a non-trivial conformal field theory. This particular property of 
non-Abelian field theories is of more than idiosyncratic interest, due to the appearance of confor- 
mality in theories beyond the standard model, such as walking technicolour [3,4]. 

1.1 Previous results 

The beta function can directly provide us with upper limit of this region, which is the per- 
turbative point where asymptotic freedom is lost. The lower limit depends on non-perturbative 
dynamics. Several analytic approaches have been applied to this problem, including renormalisa- 
tion group techniques to estimate the critical point [5] and degree of freedom analysis [6], arriving 
at Nj- < AN-\/ 1 — \6/^\N^ for an SU {N) theory with Nf flavours. Renormalisation group flow 
equations predict N'j: = 10(1) [7]. A supersymmetry inspired all-orders beta function leads to a 
bound of the conformal window [8] of A'y > 8.25 for an SU(3) gauge group. A recent study at 
finite temperature, on the base of a truncated renormalisation group flow calculation [9], found a 
value of A'^^ close to 12. 

On the lattice side, the authors of [1 1] performed a pioneering study at Nf = 10 and Nf = 12, 
concluding the former to be outside of the conformal window. Subsequent work by the Columbia 
collaboration [12] was not completely conclusive, but tentatively set Nf = 8 outside of the con- 
formal window. Iwasaki and collaborators, however, found a value as low as A'y = 6 [13]. In the 
already mentioned paper [10], the authors consider A^/ = 16, the largest value oiNf still compatible 
with asymptotic freedom, and confirm that chiral symmetry is not broken at weak coupling. Sev- 
eral other studies have observed phase transitions as a function of the gauge coupling for Nf = % 
[14, 15, 16, 17], leaving the question of the fate of these transitions in the continuum limit open. 
Recent work by Appelquist and collaborators [18], using the Schrodinger functional and a step 
scaling method, concludes that Ay should be close to 12. A new study by Fodor et al. [19], fo- 
cussing on the spectrum of low eigenvalues for 8 and 12 flavours, concludes cautiously that both 
are in the regular phase of QCD, allowing for taste breaking effects. Recently, there has also been 
active interest in the properties of theories using fermions in different representations [20, 21]. 

In summary, older studies at A^/ = 8 suggested that the theory is already in the conformal 
window, while others confirmed analytic calculations [2, 5, 8, 9], indicating that Nj- > 8, close to 
the upper hmit for A^^^ calculated in [6]. 

1.2 Our approach 

Since it is a defining property of the conformal phase, the most straightforward approach to 
determining if Nf = 8 is in the conformal window would seem to be measuring the amount of 
spontaneous chiral symmetry breaking at zero temperature. However, the strong coupling limit of 
lattice QCD is always confining, regardless of the value of Ay (see again e.g. [10]). So even for 
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those Nf in the conformal window, a phase of broken chiral symmetry exists. However, this strong 
coupUng phase is a lattice artifact and will not be traversed by any physical RG trajectories. When 
interpreting the simulation results, one therefore has to be certain that the chosen parameters did 
not put one in an unphysical phase. On the other hand, our finite lattice always introduces some 
finite temperature. Observing a chirally symmetric vacuum, therefore, could simply mean we 
are looking at a thermodynamically generated deconfined phase of what is otherwise a confining 
theory. A single observation of the chiral properties, therefore, does not actually imply anything 
yet. If a stratagem such as the one proposed here is to work, one needs to determine the nature of 
the transition itself. What we present here is a determination of the conformal qualities of eight 
flavour QCD from the nature of the observed phase transition as published previously [22]. 

2. Simulation setup 

When Nf is chosen outside of the conformal window, no transition is expected at zero temper- 
ature, i.e. when the temporal extent of the lattice is effectively infinite. If a lattice transition occurs, 
it should smoothly move towards {3 —k>o as Nt ^ °°. When it is chosen inside the conformal win- 
dow, the theory should become insensitive to infrared length scales. Any phase transition at low 
beta would have to be bulk as a consequence; it should survive the A^, — > oo limit and be insensi- 
tive to its exact value. In practise, one could distinguish between these behaviours by varying the 
values of Nt and j3 for a given theory. If the critical value of j3 moves appropriately with Nt, we 
are observing a thermodynamical transition at a well defined physical temperature = ^^jj^- If 
it remains stationary, the transition should be bulk. The former option tells us immediately that the 
theory is not conformal, the latter leaves open the possibility of a transition at weaker coupling. 

For our systematic study of the phase transition, we used two values of Nt, namely 6 and 12. A 
third calculation was done at the even smaller value of A', = 4, a region that should be problematic 
for scaling analysis due to the coarse lattice spacing at which the transition should occur there. 
The value found there was used to generate confidence in our treatment of the N = 6 results. To 
check for the influence of finite volume effects and deduce the scaling properties of observables, 
simulations were run for three spatial extents of the lattice Ns = 12,20,24 for N = 6, and A^^ = 24 
for A^, = 12. We chose a fixed value of the lattice degenerate quark mass am = 0.02 to explore the 
critical region. 

We use an improved Kogut-Susskind fermion action, the "Asqtad" action which removes lat- 
tice artifacts up to 0{a^g^), the details of which can be found in [22]. The tadpole parameter 
Mo is the only parameter in need of tuning in this action, and it can determined using its role in 
cancelUng the the Hartree term in a perturbative expansion [23]. As such, it should reflect the 
expectation value of the physical field strength. Since the improvement feeds back into plaquette 
value stably, a self consistency criterion suggests itself. This can be implemented rather efficiently 
on lattices of small sizes and interpolation with low order polynomials produces excellent initial 
guesses once several values are known. The rational hybrid monte carlo method was employed for 
the generation of ensembles. 

The chiral condensate for Nf degenerate flavours in lattice units 



Nf 



(2.1) 



4N^N 
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Figure 1: A. The chiral condensate (red) and Polyakov loop (blue) as a function of the lattice coupling )3. 
B. The scalar Xa (blue) and pseudoscalar Xk (red) chiral susceptibilities as a function of the lattice coupling 
/3, confirming the degeneracy of the chiral partners in the symmetric phase. The cumulant is also shown 
(black, top of figure). Both. Measurements were done for three different spatial extents of the lattice N^: 12 
(O)^ 20 (A) and 24 (□). The critical region has been indicated by vertical lines 

was determined by using a stochastic estimator with 20 repetitions. The chiral susceptibility, mea- 
suring the variation of the chiral condensate with varying the fermion mass % = d{\j/\i/) /dm at 
fixed j3 can be divided into a connected and disconnected component x = Zconn + ;tdisc, given in 
lattice units by 

Nf 



2 

^ ^disc 



Nj 



16N^N, 



(Tr[M-^]^)-(Tr [M 
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(2.2) 



respectively, all written in terms of traces of the staggered fermion matrix M. The connected 
and disconnected contributions to the chiral susceptibility were measured separately, in the same 
manner as [24]. 

We can use the chiral susceptibility and the chiral condensate to define two physically relevant 
quantities 



Xa = X 



1m 



and 



m 



(2.3) 



(2.4) 



They are related through Ward identities to the spacetime volume integral of the scalar (a) and 
pseudoscalar (tt) propagators and should therefore become degenerate when chiral symmetry is 
restored. Their associated cumulant = XalXn therefore a most useful physical observable in 
analysing chiral phase transitions. 

3. Results 



Results for the chiral condensate (and the Polyakov loop) for Nt = 6 lattice are shown in 
figure lA. From the small differences found at different spatial extents, it was concluded that our 
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Figure 2: The chiral condensate at A^, = 12 and Ns = 24 in lattice units as a function of the lattice coupling 
/3 . Best fit curves are superimposed and the vertical lines indicate the critical region. The absolute value of 
the finite difference between measured values of the condensate, as an approximation of its first derivative, 
is plotted in the bottom part of the figure with an arbitrary rescaling. It shows a peak at j3 = 4.34. 

results can indeed be considered infinite volume estimates for j8 < 4.1 and j8 > 4.15, and the j8 
dependence is smooth. The jump between the two branches is very clear for both observables and 
suggestive of a discontinuity. A volume dependence is seen in the critical region, as would be 
expected for a first order phase transition [25] . We placed the infinite volume limit value for jSc at 
the lower end of the designated critical region, at and is bounded by jS^ = 4.1 125(125). 

The chiral susceptibilities plotted in figure IB confirm our picture. A small amount of splitting 
remains even in the symmetric phase, which is the explicit breaking induced by the mass term. The 
cumulant Rji, as defined earlier and plotted in figure IB shows most clearly the amount of deviation 
from full chiral symmetry restoration, which lies in the order of 5%. This observable also exhibits 
a clearly discontinuous jump within the critical region. 

After having identified the transition point for the coarser lattice, we performed simulations at 
Nt = 12 with spatial volume A^^ = 24 and the same lattice mass am = 0.02, the main result of which 
is shown in figure 2. We observe a jump in the chiral condensate suggestive of a discontinuity, but 
slightly distorted and smoothed by the relatively modest spatial extent. Notice also that effects of 
explicit chiral symmetry breaking due to the non zero physical mass are to be more pronounced 
here, as the bare fermion mass is now twice as big in physical units. Using a best fit value for 
the maximum of the derivative indicates an upper bound on the value for the critical coupling of 
j8, = 4.34(4). 

It is clear that the transition has shifted as a function of the temporal extent of the lattice, which 
points towards a thermal transition rather than a bulk one. To confirm this, an asymptotic scaling 
analysis is needed in order to verify that we are actually measuring a critical temperature in the 
continuum. We do so by means of the standard relation that connects the lattice cutoff A^ to the 
gauge coupling g, aAi = R{g^) with 

R{g') = {bog')-''^^'''^e-'^^''o^\ (3.1) 
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Figure 3: A. Finite difference approximation to the absolute value of the first derivative of the chiral con- 
densate as a function of T /Tc, determined using perturbative RG scaling. Data at A^; — 6, A^., — 24, 12 (blue) 
are compared with data at A^, — 12, A^j = 24 (red). The point T /Tc — 1 corresponds to /3( (A^, = 12) = 4.34. 
A Gaussian fit is superimposed to the A^, = 12 data, while indicative Gaussian curves are shown for A', = 6. 
The baseline is subtracted and we only indicate data near the critical points. The lattice determined critical 
regions at A'r = 6 (blue) and A^, = 12 (red) are indicated by vertical lines. B. Scaling of the critical temper- 
ature presented in terms of the critical lattice coupling j3 for Ni = 6 and 12 lattices. The central line is a 
two loop asymptotically scaled line of constant physical temperature T = T^, where N, = 12, stemming from 
the lattice with lowest. All lattice results should eventually follow this line for large enough A^,, as lattice 
artifacts are increasingly suppressed by a decreasing lattice spacing. A preliminary point at A^, = 4 has been 
added for comparison. The lattice spacing here would be expected to be too course for this type of analysis, 
but the scaling relation seems to hold even up to this point, suggesting effective improvement of the action. 

where the two loop RG running of the j3 -function is accounted for, with the universal one- and 
two-loop coefficients given by 

for Nf massless flavours. Given the definition = -^^ijyj^, the scaling relation NtR{gc{Nt)) = const 
is implied. Solving for Nt = 6 and Nt = 12, we can predict ^c{gc{Nt = 6)) by knowing ^c{gc{Nt = 
12)). A strong discrepancy with the actual lattice determination might be suggestive of a bulk zero 
temperature transition, while a small discrepancy can simply be expected and imputed to violations 
of asymptotic scaling and residual effects due to a non zero fermion mass. We find an extrapolated 
value of ^^'^^ = 4.04(4), within 2% of the measured value. In terms of T /T^ this corresponds to 
good agreement for this sensitive quantity, to within 20% (see figure 3). The present results allow 
us to firmly conclude that we are seeing a true thermal transition, in other words Nf = 8 undergoes 
a chiral restoration transition at finite temperature. 

4. Conclusion 

We find pciNf = S, Nt = 6, A^, = oo) = 4.1125(125) and pc{Nf = S, Nt = 12, A^, = 24) = 
4.34(4). These values ai^e in agreement with asymptotic scaling according to the perturbative two 
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loop beta function, to about 2% in terms of J3c and to about 20% in terms of T^. This is fully 
consistent with the occurrence of a true thermal transition, confirming that SU{3) gauge theory with 
eight flavours exhibits a chirally broken phase of QCD at zero temperature and in the continuum 
limit. In agreement with recent analytical results, we conclude therefore that lower bound of the 
conformal window must lie above Nf = 8. 
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